NASA Technical Memorandum 106116 

ICOMP— 93— 12 

CMOTT-93-6 _ . 


/N-3*/ 



Remarks on Turbulent Constitutive Relations 


Tsan-Hsing Shih 

Institute for Computational Mechanics in Propulsion 
and Center for Modeling of Turbulence and Transition 

Lewis Research Center 

Cleveland, Ohio 


and 

John L. Lumley 
Cornell University 
Ithaca, New York 


(NASA— TM- 106116) REMARKS ON N93-27010 

TURBULENT CONSTITUTIVE RELATIONS 
(NASA) 12 p 

Unclas 


G3/34 0167256 


May 1993 






I j Lewis Research Center Vi \ 

'w/mrvn 


ICOMP 

i -z\ CiMWwtem 1^ 

\ q\ ReMrv« University Jr^ j 

Ohio Aerospace /J? > 
7)\ Instttuts y 

N1 S ^ 



Remarks on Turbulent Constitutive Relations 
Tsan-Hsing Shih 

Center for Turbulence Modeling and Transition (CMOTT) 
Ins titute for Computational Mechanics in Propulsion (ICOMP) 
NASA Lewis Research Center, Cleveland, Ohio 

John L. Lumley 

Cornell University, Ithaca, New York 


The paper demonstrates that the concept of turbulent constitutive relations introduced 
by Lumley (1970) can be used to construct general models for various turbulent corre- 
lations. Some of Generalized Cay ley- Hamilton formulas for relating tensor products of 
higher extension to tensor products of lower extension are introduced. The combination 
of dimensional analysis and invariant theory can lead to “turbulent constitutive relations” 
(or general turbulence models) for, in principle, any turbulent correlations. As examples, 
the constitutive relations for Reynolds stresses and scalar fluxes are derived. The results 
are consistent with ones from RNG theory and two-scale DIA method, but with a more 
general form. 

I. Introduction 

Lumley 1 discussed some possible constitutive relations for turbulent correlations. One 
of his conclusions was that turbulent constitutive relations may exist under situations in 
which the length and time scales of turbulence are smaller than those in the mean flow field 
so that the effect of initial and boundary conditions on the turbulence is not significant. 
For example, if an unbounded mean flow evolves slowly in both time and space, there 
may exist a constitutive relation for turbulent correlations (e.g., Reynolds stresses) which 
may help to solve the turbulence closure problem in this particular situation. However, in 
most practical situations, the scales of turbulence are of the same order of magnitude as 
those of the mean flow field, hence so called turbulent constitutive relations do not exist, 
and the turbulence closure problem cannot be theoretically solved. The formally derived 
“turbulent constitutive relations” become at most approximations (turbulence models) 
to the turbulent correlations whose validity varies from flow to flow. Regardless of the 
existence of turbulent constitutive relations, the formal procedure of deriving constitutive 
relations developed in continuum mechanics provides an useful tool for developing general 
turbulence models. 
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A constitutive relation developed in continuum mechanics is always of an equilibrium 
and local form, that is, it is a function of arguments at the present time and the local 
point. This is, of course, not true for most turbulent flows. For a general turbulent flow, 
a constitutive relation should contain the time history and spatial information for the 
arguments in question. However, as a model approximation, we always neglect time and 
spatial effects and consider the relationship at the present time and the local point as the 
first order approximation in the time and spatial expansions of the functional form. 


The procedure for deriving a constitutive relation usually includes two steps. The first 
step is to make sin assumption about the general relationship between the turbulent corre- 
lation and other “known” quantities (which are the arguments of the function) using the 
information from whatever observations or experience we have. At this stage, dimensional 
analysis (t theorem) can be used to regroup and reduce the number of independent argu- 
ments in question. The second step is to determine the detailed form of the relationship 
using invariant theory. 

In the application of invariant theory, the most important procedure is to determine the 
number of independent invariants. This determination will become more complicated 
if the arguments in question have more than two tensors. The number of independent 
invariants depends on the number of independent tensors which it is possible to form using 
the arguments in question. In this paper, we shall show a procedure for determining the 
independent tensors which can be formed by two general tensors using generalized Cayley- 
Hamilton formulas. Thereafter, as an example, we shall briefly derive a constitutive relation 
for the turbulent Reynolds stresses. A possible constitutive relation for the scalar flux is 
also given. This approach can be extended to more complicated turbulent correlations, for 
example, the high order turbulent correlations appearing in the transport equations of the 
second order moments. 


II. Generalized Cayley-Hamilton formulas 

Rivlin 2 showed that there are several generalized Cayley-Hamilton formulas relating ma- 
trices (product of several matrices A,B,C ...) of higher extension to matrices of lower 
extension. Here, we follow Rivlin to derive some relations of this type which will be used 
latter. 
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One of Rivlin’s basic relations is 

(ABC + ACB + BCA + BAC + CAB + CBA ) - A(trBC - trB trC ) 

- B(trCA - trCtrA ) - C(trAB - trA trB) 

- ( BC + CB)trA - (CA + AC)trB - (AB + BA)trC (1) 

— I (trA trB trC — trA trBC — trB trC A 

-trCtrAB + trABC + trCBA) = 0 

where I is the unit matrix or identity tensor and tr represents the trace operator. 

Replacing C with A in Eq.(l), we obtain 

ABA = -A 2 B - BA 2 + A(trAB - trA trB) 

+ - B(trA 2 - trA trA) + (AB + BA)trA + A 2 trB (2) 

2 

+ I[trA 2 B - trA trAB + trB (trA trA - trA 2 )} 

Replacing C with B in Eq.(l), we obtain 

BAB = -B 2 A - AB 2 + B(trBA - trB trA) 

+ i A(trB 2 — trB trB) + (BA + AB)trB + B 2 trA ^3) 

2 

+ I[trB 2 A - trB trB A + \trA(trB trB - trB 2 )] 

Zt 

Eqs.(2) and (3) indicate that the matrices ABA and BAB of extension 3 can be expressed 
by polynomials of matrices of extension 2 or less. 

Multiplying Eq.(2) by A both on the left and on the right, and then adding the results, 
we obtain 

ABA 2 + A 2 BA =ABA trA + A 2 trAB + A(trA 2 B — trA trAB) 

— B detA + I detA trB 

where det is the determinant operator. Multiplying Eq.(3) by B on the left and 
right, and then adding the results, we obtain 

BAB 2 +B 2 AB =BAB trB + B 2 trB A + B(trB 2 A - trB trB A) 

— A detB + I detB trA 

It follows from Eqs.(2) and (3) that the right hand side of Eqs.(4) and (5) is a polynomial 
of matrices of extension 2 or less. 


( 4 ) 

on the 

( 5 ) 
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( 6 ) 


Replacing B with B 2 in Eq.(4), we obtain 

AB 2 A 2 + A 2 B 2 A =AB 2 A It A + A 2 trAB 2 + A{trA 2 B 2 - trA trAB 2 ) 

- B 2 detA + I detA trB 2 

Replacing A with A 2 in Eq.(5), we obtain 

BA 2 B 2 + B 2 A 2 B =BA 2 B trB + B 2 trB A 2 + B{trB 2 A 2 - trB trB A 2 ) 

— A 2 detB + I detB trA 2 

Replacing B with B 2 in Eq.(2), we obtain 

AB 2 A = —A 2 B 2 - B 2 A 2 + A(trAB 2 — trA trB 2 ) 

+ - B 2 {trA 2 - trA trA) + {AB 2 + B 2 A)trA + A 2 trB 2 (g) 

2 

+ I[trA 2 B 2 - trA trAB 2 + ]-trB 2 (trA trA - trA 2 )) 

& 

Replacing A with A 2 in Eq.(3), we obtain 

BA 2 B = -B 2 A 2 - A 2 B 2 + B(trBA 2 - trB trA 2 ) 

-f — A 2 {trB 2 - trB trB) + {BA 2 + A 2 B)trB + B 2 trA 2 (9) 

2 

+ I[trB 2 A 2 - trB trB A 2 + \trA 2 {trB trB - trB 2 )) 

Eqs.(8) and (9) indicate that the matrices AB 2 A and BA 2 B of extension 3 can be expressed 
by polynomials of matrices of extension 2 or less. As a result, the right hand sides of Eqs.(6) 
and (7) are also polynomials of matrices of extension 2 or less. 


III. Number of independent tensors formed by two general tensors 

Now let us show that the number of independent tensors formed by two general tensors A 
and B is 18. 

Rivlin 2 showed that any matrix product in two 3x3 matrices may be expressed as a 
polynomial of these matrices of extension 4 or less. For example, suppose we have a 
matrix product II of extension 5: 


n = aba 2 b 2 a 


( 10 ) 
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This can be written as 


n = AC A 


( 11 ) 


where C = BA 2 B 2 . From Eq.(2), II may be viewed as a polynomial of matrices in A and 
C of extension 2 or less. Note that C itself is a matrix in A and B of extension 3, therefore, 
II may be expressed by a polynomial in A , B of extension 4 or less. Therefore, we only 
need to consider the possible matrices (or tensors) of extension 4 or less formed by A and 
B. 


We now show that there are only two independent tensors of extension 4. The possible 
tensors of extension 4 are the following 8 tensors: 


ABA 2 B 2 , BAB 2 A 2 , A 2 BAB 2 , B 2 ABA 2 , 
AB 2 A 2 B, BA 2 B 2 A, A 2 B 2 AB, B 2 A 2 BA. 


With Eqs.(4) and (2), A 2 BAB 2 can be expressed by ABA 2 B 2 plus a polynomial in A and 
B of extension 3 or less. Similarly, with Eqs.(5) and (3), B 2 ABA? = —BAB 2 A 2 + poly- 
nomial in A and B of extension 3 or less, with Eqs.(7) and (9), AB 2 A 2 B = ABA B ■+■ 
polynomial in A and B of extension 3 or less, with Eqs.(6) and (8), BA 2 B 2 A = —BAB 2 A 2 
+ polynomial in A and B of extension 3 or less. From Eqs.(4) and (5), we may show that 
jp B 2 AB = ABA 2 B 2 + polynomial in A and B of extension 3 or less, and that B 2 A 2 BA 
= BAB 2 A 2 + polynomial in A and B of extension 3 or less. Therefore, only two tensors 
of extension 4 in Eq.(12) are independent, we select them as 

ABA 2 B 2 , BAB 2 A 2 (13) 


Now we show that there are only four independent tensors of extension 3. The possible 
tensors of extension 3 axe the following 8 tensors: 

ABA 2 ,A 2 BA, BAB 2 , B 2 AB, AB 2 A 2 , A 2 B 2 A, BA 2 B 2 ,B 2 A 2 B (14) 

Using Eqs.(4), (5), (6) and (7), we find that only four of them in Eq.(14) are independent. 
Let us select them as 

ABA 2 ,BAB 2 ,AB 2 A 2 ,BA 2 B 2 (15) 


There are eight possible tensors of extension 2 which are all independent: 

AB, BA, AB 2 , B 2 A, 

A 2 B, BA 2 , A 2 B 2 , B 2 A 2 . 


(16) 
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( 17 ) 


There are four possible tensors of extension 1 which are also independent: 

A, A 2 , B , B 2 . 

Therefore, we have proved that only 18 tensors can be formed independently by two general 
tensors. 


IV. Constitutive relation for Reynolds stresses UiUj 


We assume 

UiUj = FijiUij, k , f) 

which indicates that the turbulent stresses depend on the mean velocity gradient Uij and 
the scales of turbulence characterized by the turbulent kinetic energy k and its dissipation 
rate e. Applying the 7r theorem, the arguments can be regrouped as 





Accordingly, we may write 


UiUj 

2k 


= FijiAij) 


Noting that Aij is a general tensor, if we define its transpose as 


B„ = A% = (18) 

then, B ^ A. In order to obtain a general relationship, we first look for a general tensorial 
form of Fij(Aij,Bij), then replace Bij by Af-. We have proven that the independent 
tensors formed by A and B are the following 18 tensors: 


a,a 2 ,b,b 2 , 

AB , BA, AB 2 ,A 2 B, BA 2 ,B 2 A, A 2 B 2 ,B 2 A 2 , 
ABA 2 , BAB 2 , AB 2 A 2 ,BA 2 B 2 , 
ABA 2 B 2 ,BAB 2 A 2 . 


Following Lumley 3 and applying invariant theory, we may obtain 

— = a\Sij + a^A + a^B + U4-A 2 + a$B 2 

2k 

+ clqAB + a?BA + a&AB 2 + a$A 2 B 
•f ujo BA 2 oi^B 2 A o\%A? B 2 “l" u\$B A 
+ a\\AB A? + a,i5 B AB 2 + ui$AB 2 A 2 + anBA 2 B 2 
+ a ls ABA 2 B 2 + ai*BAB 2 A 2 
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where aj — a* 9 are scalar functions of the invariants of the tensors in question. 


Using the conditions: UiUj = UjU,, UiUi — 2 k and Bij = Afj = Aji , we obtain 


02 = 03, 04 = 05, a 8 = o 9, oio = on, 

Ol4 = O15 = 016 = Oj7 = Oj8 = Oj9 = 0 


( 21 ) 


and 


o 1 


Therefore, we obtain 


— — (l — 2o2^tt — 2a^AijBij (oj *1" Ori^AijBji 
O 

- 2 (a 8 + a 10 )AijBji — (012 + a^A^B^] 


u iUj = + 2a2—^-(Uij + Uj t i — -U iAtj ) 

+ 2 a 4 ^-(Ui,kUk,j + Uj,kUk,i - 

+ 2a t ~(V i , k U j , k - |fi«y) 

+ - ifiiy) 

+ 2a.£-(Vi, k Ul t + Ufj,V i, k - |fi«y) 

+ 2 a^{U k ,iUl,i + UkjUli - 
+ 2 a l 2 ^-(ul k ul„-^ns ii ) 

+ 2 a lz ^-(Ul, i Ul,i-l^ii) 


(22) 


( 23 ) 


where 

n = u iJ u jli , n = = n = (24) 

It is noticed that the first two terms on the right hand side of Eq.( 23 ) represent the 
standard k-e eddy viscosity model, and that the first five terms of Eq.( 23 ) are of the same 
form as the models derived from both the two-scale DIA approach (Yoshizawa 4 ) and the 
RNG method (Rubinstein and Barton 5 ). 

Eq.( 23 ) is a general model for uiuj. It contains 8 undetermined coefficients which are, in 
general, scalar functions of various invariants of the tensors in question, for example, SijSij 
and toijflij which are (fi + H)/2 and (fi - II )/2 respectively. The detailed forms of these 
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scalar functions must be determined by other model constraints, for example, realizability, 
and by experimental data. Eq.(23) contains 9 terms; however, its quadratic tensorial form 
may be sufficient for practical applications, expecially when \Ui t j\k/ e is less than one. 


V. Constitutive relation for Scalar Flux 9u{ 


We assume 

hr i =F i {U iti ,T,i,k i e,0\e e ) (25) 

where & 1 is the variance of a fluctuating scalar and e$ is its dissipation rate. Eq.(25) 
indicates that the scalar flux depends on not only the mean scalar gradient Tj, but also 
the mean velocity gradient Ui t j and the scales of both velocity and scalar fluctuations 
characterized by k, e, 0 2 , e$. 


Applying the i r theorem, we have: 


Mi- e, t - ( ££<? ) 1/3 r .i* 


k , 0 2 
r = - — 
£ £0 


Accordingly, we may write 


J Ui , = Fi{Aij ,0,i, r) 

( 02fc )l/2 

Using arg um ents similar to the ones made above for Reynolds stresses, we obtain 


+ + aiUj.kUk.i + + ajUkjUkj ) Tj 

£ 2 £ £0 

+ ) 1 / 2 (a 8 Ui t kU 2 k + a 9 U 2 k Uj,k + aio Uk,iU k j + a u Ul ti Uk,j ) Tj 

£ a £ £0 

& t 

+ i,kV e,kUfj + ois Uk,iUk,iUj tt ) Tj 

+ iUvVh.Vh + avUkM/H' 


N £ e$ 


+ ai 6 Ui t kUt t kUj <rn U 2 tm + a 19 U k ,iUk, t U 2 m>i U 2 m>j ) Tj 


(26) 


(27) 


The coefficients a\ - aig are, in general, functions of the time scale ratio r and the other 
invariants of tensors in question. 
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If we assume the time scale ratio r « constant, for example, 0 2 /ee « k/e, then Eq.(27) 
becomes 


9ui = <*i — T ti + —(<*2 Uij + a 3 Uj,i)Tj 

£ £ 

+ Ui t kUk,j + asUj,kUk,i + a’6Ui t kUj t k + aiUk,iUk,j ) Tj 


+ 

+ 

+ 


■^{ a sUi,kUj k + asUfkUj'k + autU^iUlj + <*n Ul ti Uk,j ) T t j 

^5 (<*12 Uf tk Ulk + <*1 zUl yi U kt j 
+ auUi'kUt'kUij + a 15 Uk,iUk,iUj tl ) T,j 

^-(<*16 Ui,kUt,kUt,j + <*1 7Uk,iUk,t U lt 

+ CLigUi t kUi t kUi t7n Uj im + aivUkjUktiU^iU^j ) Tj 


(28) 


It is interesting to note that the conventional eddy viscosity model for the scalar flux is 
just the first term on the right hand side of Eq.(27) or (28), and that the models derived 
from two-scale DIA (Yoshizawa®) and RNG method (Rubinstein and Barton"^) are the first 
two terms of Eq.(28). 


VI. Conclusion 

We have demonstrated that the combination of dimensional analysis and invariant the- 
ory is a powerful tool for developing “turbulent constitutive relations” (or general turbu- 
lence models) for various turbulence correlations. The way of forming generalized Cayley- 
Hamilton for mulas for determining independent tensors is also shown. As examples, the 
general turbulence models for Reynolds stresses and scalar fluxes are derived. The results 
from RNG theory and the two-scale DIA method are the first few terms of the constitu- 
tive relations derived in this paper. This technique can be extended to other turbulence 
correlations including higher order correlations appearing in second order closures. 
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